Effect of quantum interference in the nonlinear conductance of microconstrictions 
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The influence of the interference of electron waves, which were scattered by single impurities, on 
nonlinear|-(juantum conductance of metallic microconstrictions (as was recently investigated experi- 
mentally Lfu ) is studied theoretically. The dependence of the interference pattern in the conductance 
G (V) on the contact diameter and the spatial distribution of impurities is analyzed. It is shown 
that the amplitude of conductance oscillation is strongly depended on the position of impurities 
inside the constriction. 
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I. INTRODUCTION. 



The Scanning Tunneling Microscopy (STM) and the 
Mechanically Controllable Break-junction (MCB) tech- 
niques offer an opportunity to study the conductance of 
metallic contacts consisting of only a few atoms (quan- 
tum contacts). The electrical conductance of such con- 
tacts, at small bias voltage is proportional to the number 
of propagating electron modeg., TV, each one contributing 
an amount of Go = 2e^ /h U. With increasing the di- 
ameter of the contact the energies of modes continuously 
decreases, but the number of modes increases whenever 
a new mode fits into the constriction cross section. This 
number, N, is limited by the requirement that the ki- 
netic energy for the transverse motion is smaller than the 
Fermi energy ep- When a new mode is occupied, a new 
quantum channel is opened. The conductance then un- 
dergoes a jump of Gq. Such quantization of conductance 
has been observed in both two and three dimensional con- 
tacts with diameters comparable to Fermi|-.wave-length 
\f — h/pp {pp is the Fermi momentum)a~0 . Jumps in 
conductance are also expected to occur, at the constant 
contact diameter while bias voltage is varied. If the bias 
eV is larger than the distances between the energy levels 
of quantum modes, it is possible to change the number 
of opened modes by changing the voltage V. At a cer- 
tain threshold voltage a channel is opened or closed for 
one direction of the electron wave vector along the con- 
striction and consequently conductance suffers a Go/2 
stepwise changeO. 

Quantum interference effects have been studied in dif- 
ferent mesoscopic systemalJ. In ballistic metallic micro- 
constrictions it manifests itself as fluctuations in conduc- 
tance \dien a magnetic fieldtll or an electrical voltage is 
appliedEj. Now experimental efforts have been done us- 
ing MCB techniques, in order to measure conductance 
as a function of voltage in atomic-size point contactso. 
A prominent feature of these measurements, is the exis- 
tence of small random voltage-dependent fluctuations in 



conductance, far from steps. The measurements§ clearly 
indicate suppression of the fluctuations for conductance 
values near the integer multiples of the conductance 
quantum. Similar results have been reported by using a 
STM to show the strong voltage dependence of condup.- 
tance of one-atom contacts at different temperaturcffl. 
It is generally believed that the observed oscillations in 
conductance are due to the quantum interference effectsEI. 
Ludoph and co-authorsu, propose the following interpre- 
tation: The electron wave transmitted through the con- 
tact is backscattered to the contact by an impurity and 
then partially reflected at the contact. These waves in- 
terfere and change the total conductivity. The energy 
and thus the wave number of an injected electron into 
the channel, depends on the voltage. Consequently the 
interference pattern in conductance oscillates as electron 
wave number varies with the voltage. 

Although the theory developed by Ludoph et al.0 can 
explain the general feature of fluctuations but here we 
try to examine a different mechanism. Impurities (or de- 
fects) are assumed to be located inside the constriction, 
and the interference is effectively between waves scat- 
tered from the impurities. The existence of a few defects 
or impurities inside the constriction is rather natural con- 
sidering the way the contact is formed. Using the model 
of a long microconstriction we can find the conductivity 
analytically. We discuss the theory of nonlinear electron 
transport through a mesoscopic microconstriction with a 
few impurities. We show that the nonlinear dependence 
of the quantum conductance on the voltage is obtained 
from this model. The form of this dependence is affected 
not only by the distances between impurities, but also by 
their positions inside the constriction. 

In Sec. II the model Hamiltonian is discussed and is 
used to obtain a general expression for the nonlinear con- 
ductance. In Sec. Ill a ^-function potential is assumed 
for the interaction of electrons with impurities and a sim- 
plified equation for the conductance is obtained. Within 
the framework of perturbation theory, a general analyt- 
ical equation is also derived for conductance of the sys- 
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tern, for arbitrary number of quantum modes and arbi- 
trary number of impurities located in arbitrary positions. 
These analytical results are illustrated by numerical cal- 
culations for the contact in the form of a long cylindrical 
contact. A brief discussion of result is given in Sec. IV. 



II. GENERAL EQUATION FOR THE 
NONLINEAR CONDUCTANCE OF THE LONG 
QUANTUM MICROCONSTRICTION. 



Let us consider a long narrow constriction, which con- 
nects two bulk metals, assuming eV <C ep- The ge- 
ometry is shown in Fig. 1. We assume that the con- 
tact shape is smooth on the scale of the wavelength Xp- 
This condition assures that different transverse modes 
pass through the ballistic contact independently (adia- 
batic approximationllj) . We also assume that the con- 
tact length is much larger than its diameter and we can 
neglect the constriction end effects. Under these approxi- 
mations, the electrical field inside the contact far from the 
ends is negligible and the energy e of ballistic electrons 
deppads. only on the sign of velocity along the contact 



axis 



The Hamiltonian H of the electrons contains 



the following terms: 



H — Hq + Hi + Hi, 



where 



is Hamiltonian of free electrons, and 
n \- ■ t 



(1) 



(2) 



(3) 



describes the influence of applied bias voltage. Hint de- 
notes interaction of electrons with impurities, and de- 
pends on the positions of impurities in the constric- 
tion; the operator (cq) creates (annihilates) a conduc- 
tion electron with the wave function (/Jq,, and energy e^. 
The electron wave functions and eigenvalues are 



(4) 
(5) 



Lpair) = V'/3(R)exp ( -p,_z 



2me 



where a = (/3,pz) , /? is the set of two transverse quan- 
tum numbers; is the momentum of an electron along 
the contact axis, r = (Ti,z) , R is a coordinate in the 
plain, perpendicular to the z axis me is the electron mass. 

In zero approximation in Hint the current Jq through 
the contact area Sc is 



Jo = eScTr (v^po) 



(6) 



where 



Po = fF{Ho + m); 



(7) 



Vz = Pz/'Tn is the electron velocity; fp is the Fermi func- 
tion. Using the Eqs.(6),{7), and wavefunctions (4), we 
find the equation for the ballistic conductance: 



G 



1 - H 



eV 



(8) 



At zero temperature and y — > 0, this formula de- 
scribes the well known Go steps of quantum conductance 
and in the quasjclassical case it turns into the Sharvin 
conductanceEErM. 

In order to investigate the influence of single impuri- 
ties on the nonlinear quantum conductance of the point 
contact, we ijse the method, which was developed Kulik 
and otheralflo. The change in the electrical current A J 
is related to the rate of energy dissipation by the relation: 



dE _d {Hi) 
'dt ~ dt 



(9) 



Differential of (Hi) with respect to time t is we ob- 
tained from Heisenberg equation. The change A J of the 
current due to interactions of electrons with impurities; 
would then be 



where 



AJV ^^{[Hi{t),H,nt (t)]), 
in 



(...)=Tr(pW...). 



(10) 
(11) 



All operators are in the interaction representation. 
The statistical operator p (t) satisfles the equation 

dp 



ifi- 



dt 



[Hrnt (t) , p (t)] , 



(12) 



which can be solved using perturbation theory in Hint 
(but for the arbitrary Hi): 



P (t) = Pq + 



^ J dt' [H,nt it'),Po] 



(13) 



We would then have 
t 



AJ = - 



1 



dt'Tr{po[[Hi,H,nt{t)],H,ntit')]). 



The decrease in total conductance AG ^ 
the quantum interference is deflned as 

dAJ 



Go = — 



dV 



(14) 

-G2, results 
(15) 



If the applied bias eV is much smaller than the differences 
between the energies of modes, the Eq. (15) describes 
the dependence of total conductance on the voltage V. 
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III. THE CONDUCTANCE OSCILLATIONS. 



Now using the general Eqs. (14), (15), we investigate 
the behaviour of G2 for the case of 5— function scattering 
potential. The Hamiltonian Hint can be written as 



H^ntijj) = g^ IYj) V3a' (r^) ^C^' 



(16) 



Here g is the coupling constant of the interaction of an 
electron with an impurity located in the point Vj . 



cos 



cos 



Substituting Eqs. (7), (16) into Eq.(14), after some sim- 
ple but cumbersome calculations we find 



A J = --^9^ X! X! i'^Wnvza - signvza') x 

a, a' i.j 

</'a'(rj)<^a(r.)<^Q'(ri)<^a(l"j)(/a' " - £«), (17) 



where = jp{e^ ^signVza) ■ At zero temperature 
fF = Q {ep — s) , the Eq.(17) can be further simplified. 
Using the wave functions (4), we obtain for nonlinear 
part of conductance the following equation: 



(+) 



Az 



WW 



6 (ep - Eft 



(18) 



1-) 



e {SF - Efi,,, - ^)} 7^;,(R,)V';(R^)#'(R.)V'^(RJ)• 



Where 



e/3 



(19) 



The cosine terms in the Eq. (18) describe the con- 
ductance oscillations due to the interference of electrons 
waves scattered by impurities. The transverse parts 
'0/3 (Rj) of wave functions contain the mesoscopic effect 
of impurity positions inside the constriction. The equa- 
tion (18) diverges at p^^^ — 0. Physically it means that 
in the Born approximation the slowly moving electron is 
repeatedly scattered on the impurity. In this case the per- 
turbation theory (Born approximation) is not valid any 
more, and we must take into account the interference of 
partial waves under the electron scattering by impurity. 
We assume that energy levels are not very close to the 
boundary energies sp ± ^ and the quantity G2 added 
to the ballistic conductance Gi (8) is small. 

For the numerical calculations we have used the free 
electron model of point contact consisting of two infinite 
half-spaces connected by a long ballistic cylinder of a ra- 
dius R and a length L (Fig.l). 

In the limit L 00 the electron wave functions ipa' (r) 
and energies Eq, can be written as 



ifia, (r) 



m+1 



exp [ imip + j^PzZ 
(20) 



channel, 7mn are the n-th zero of the Bessel function J^. 
Since the electron energy has degeneracy for azimuthal 
quantum number m ( as a result of the symmetry of the 
model), quantum modes with ±m have the same contri- 
bution in conductance. In this model, Sharvin ce«duc- 
tance has not only steps C?Oj but also steps 26*0113. In 
Fig. 2 the dependence of the interference pattern on the 
number of impurities inside a constriction with constant 
radius is shown. It shows that as a result of the inter- 
ference of electron waves, which were scattered by differ- 
ent impurities, the interference maxima in G2 {V) depen- 
dence, may both be depressed and increased. The inter- 
ference oscillation of the conductance depends strongly 
on the number of opened quantum modes which follows 
from the dependence of its maximum value of longitudi- 
nal electron momentum (see, Eq. 3) on the contact size. 
The voltage dependence of G2 for different contact sizes 
are shown in Fig. 3. Fig. 4. illustrates how the chang- 
ing in the nonlinear dependence G2 {V) changes with the 
contact size. It corresponds to the case of two impurities 
in the contact. The position of impurities and the num- 
ber of opened quantum modes are kept constant. The 
difference in the interference oscillations is a result of the 
changing in the relative positions of nodes and maxima 
of the electron wave function from the points, in which 
the impurities are situated. 



IV. CONCLUSION 



where 



2mp 



2mei?2 



7' 

Imn 



(21) 



We have used cylindrical coordinates r — (p, ip, z) with 
z along the axis of cylinder. Here a = {n,m,pz) are 
the quantum numbers, fl = ttR^L is the volume of the 



The dependence of quantum conductance of metallic 
ultrasmall contacts containing impurities on bias voltage 
has been theoretically studied. We have shown that im- 
purities situated inside the quantum microconstriction 
produce a nonlinear dependence of the conductance on 
the applied voltage, which is the result of the interference 
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of electron waves reflected by impurities. The transmis- 
sion probability of the electron through constriction de- 
pends on the relation between the electron wave length 
A and Azy, the projection of distances between impuri- 
ties along the channel. It is maximum when the condition 
Azy = ^ ( n is integer) is satisfied. Since the electron 
momentum depends on the applied bias, one can change 
the transmission by changing the voltage. Our numer- 
ical calculations show that the resulting nonmonotonic 
dependence of the conductance, js, similar in shape to 
the ones observed in experimentldfl. The amplitude of 
interference pattern is sensitive to the transversal posi- 
tion of impurities inside the constriction. If the impurity 
is located near the point, where the electron wave func- 
tion corresponding to the nth quantum mode vanishes, 
then the decreasing of transmission of that mode would 
be negligible. 

We acknowledge fruitful discussion with M.R.H. Kha- 
jehpour and I.K.Yanson. 



E.N.Bogachek, A.M. Zagoskin, I.O. Kulik, Sov. J. Low 

Temp. Phys. 16, 796 (1990). 

I.O. Kulik, JETP, Pis'ma, 5, 423 (1967). 

I.O. Kulik, A.N. Omelyanchouk, I.G.Tuluzov, Sov. J. Low 

Temp. Phys., 14, 149 (1988). 

Yu. V. Sharvin, Sov. Phys. JETP, 21, 655 (1965). 



^ B.Ludoph, M.H. Devoret, D.Esteve, C.Urbina and J.M. 

van Ruitenbeek, Phys. Rev. Lett. 82, 1530 (1999). 
^ C. Untidt, G. Rubo Bollinger, S. Vieira, N. Agrait,In Qun- 

tum and Mesoscopic Phenomena and Mesoscopic Devices 

tn Mtcroelectromcs (Bilkent University, Ankara, Turkey, 

1999) pp 78. 

^ For review, see Solid State Physics 44, Edit. H. Ehrenreich 

and TurnbuU (Academic Press, New York, 1991) pp 1-228. 
* B.J. van Wees, H. van Houten, C.W.J. Beenakker, J.G. 

Williamson, L.P. Kouwenhoven, D. van der Marel and C.T. 

Foxon, Phys. Rev. Lett. 60, 848 (1988). 
^ D.A. Wharam, M. Pepper, H. Ahmed, J.E.F. Frost, D.G. 

Hasko, D.G. Peacock, D.A. Ritchie, G.A.G. Jones, J. Phys. 

C: Solid State Phys. 21, L209 (1988). 
® J.I. Pascual, J. Mendez, J.Gomez-Herrero, A.M. Baro, 

N.Garsia, Vu Thien Binh, Phys. Rev. Lett. 71, 1852 (1993). 

N. Agrait, J.G. Rodrigo, S. Viera, Phys. Rev. B 47, 12345 

(1993). 

® C.J. MuUer, J.M. van Ruitenbeek, L.J. de Jongh, Phys. 

Rev. Lett. 69, 140 (1992). 
^ A.M. Zagoskin, JETP Lett., 52, 435 (1990). 

Y. Imry, Introduction to Mesoscopic Physics (Oxford, New 

York, 1997). 

" P.A. Lee, A.D. Stone, Phys. Rev. Lett. 55, 1622 (1985). 

P.A.M. Holweg, J. A. Kokkedee, J. Caro, A.H. Verbruggen, 

S. Radelaar, A. CM. Jansen, and P.Wyder, Phys. Rev. 

Lett. 67, 2549 (1991). 
" L.I. Glazman, G.B. Lesovik, D.E. Khmel'nitskii, R.I. 

Shekhter, JETP Lett. 48, 238 (1988). 
" I.O. Kulik, A.N. Omelyanchouk, R.I. Shekhter, Sov. J. Low 

Temp. Phys. 7, 740 (1977). 
^^l.O. Kulik, I.K. Yanson, Sov. J. Low Temp. Phys., 4,596 

(1978). 

A. Szafer, A.D. Stone, Phys. Rev. Lett. 62, 300 (1989). 



4 




I I I I I I I I I I I I I I I I I I I I I I 1 

0.05 0.1 0.15 0.2 eV 

(Fermi energy units) 



FIG. 2. The normalized conductance G2/G1 as a function of voltage in a three modes channel {2'jtR = 4:.2Xp,) for 
the different number k of impurities; A; = 2 for solid line, fc = 3 for dashed line and fc = 4 for dotted line. 
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FIG. 3. The normalized conductance G^jGx as a function of voltage for a channel with two impurities at different 
number of opened quantum modes; single mode {^-kR = 2.7 Xp) , for solid line, three modes {2itR = 4.2Air) , for dashed 
line and five modes {2'jtR = 5.5Air) for dotted line. 




FIG. 4. The changing of the interference pattern in the G2 {V) dependence of three mode channel, which contains 
two impurities, with increasing the contact diameter. The distance between impurities and its distances from the 
contact axis is the same for all values of R. The radius R is changed in the interval, in which the number of opened 
quantum modes is constant. 
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